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' A perturbative study of the Holstein Molecular Crystal Model which accounts for lattice structure 

I and dimensionality effects is presented. Antiadiabatic conditions peculiar of narrow band materials 

and an intermediate to strong electron-phonon coupling are assumed. The polaron effective mass 
I depends crucially in all dimensions on the intermolecular coupling strengths which also affect the 

size of the lattice deformation associated with the small polaron formation. 
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I.INTRODUCTION 



O i' Several studies published over the last years have ad- 
^ [ dressed the questions related to the existence of small 
■ polarons with itinerant properties in real systems 
I This issue is central for a possible description of high Tc 
superconductivity in terms of (bi)polaronic models 
In spite of being well defined quasiparticlcs, small po- 
'"^ larons may loose their mobility either because of a dy- 
^ ■ namical dephasing between the charge carriers and their 
O ' surrounding deformation field or because of the heavi- 
ness of the effective mass. These effects could however 
differ significantly according to the regime (adiabatic or 
. antiadiabatic) and the strength of electron- phonon cou- 
pling characterizing the system . Theoretical investi- 
, gations start generally from the Holstein Molecular Crys- 

• tal Model [^, a fundamental tool which has revealed a 
[ rich variety of behaviors in the polaron landscape through 

• the use of quantum Monte Carlo [ 12|Jl3| , density matrix 
' renormalization group techniques 14[|, variational meth- 
. ods ||7|,p^ |l7|, c luster solutions p8-|21[| and perturbative 

> ■ approaches ]22|-|2^. A numerical study of the polaron 
^ bandwidth in the first order of perturbative theory has 
y proved that the phonon momentum dependence is a key 
feature of the Holstein Hamiltonian and that the lattice 
^ ' dimensionality strongly influences the bandwidth values 
Q ] . Unlike other properties such as ground state energy 

Q : and effective mass, the bandwidth is not affected by sec- 
^ ' ond order corrections and therefore it provides a useful 
• i-H , testing bench for alternative, non perturbative attacks 
on the polaron problem [p6| . Being aware of the impor- 
^ , tance that the intermolecular forces have in the narrow- 
■ - - I ing of the polaron band, we report here on a perturbative 
numerical study of the mass enhancement in the strong 
coupling and antiadiabatic regime. The reasons why I 
choose to start from this regime are threefold: i) it is the 
easiest in the sense that the lattice deformation follows 
coherently the charge carriers and the abovemcntioned 
dephasing features can be ruled out, ii) the unit compris- 
ing electron and phonon dressing is a stable small polaron 
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that is, the size of the quasiparticle is not significantly 
broadened in some portions of our parameter space, iii) 
this regime is relevant to several classes of narrow band 
materials whose charge carriers effective masses deserve 
accurate estimates. Although the present work assumes 
that the carriers are coupled to bosonic degrees of free- 
dom having a vibrational origin, antiadiabatic conditions 
are likely to occur in excitonic systems where the charac- 
teristic frequency uj could be easily of order of leV and 
the boson can therefore follow the electron eesentially 
without retardation. In these cases the carriers effective 
mass is expected to be only moderately enhanced with 
respect to the bare electron mass. 

In Section II, the dispersive Holstein model is briefly 
reviewed while the numerical results are displayed in Sec- 
tion III. Some conclusions are drawn in Section IV. 



II. HOLSTEIN MODEL WITH DISPERSIVE 
PHONONS 



My starting point is the real space - momentum space 
representation of the Holstein Hamiltonian which reads 

i^Lj j k 

+ 5Z E + al^Jexpiik ■ r,) (0.1) 

^ k j 

c| ( Ci ) creates (destroys) a tight binding electron at the 
i — th molecular site and t is the hopping integral related 
to the bare electron half bandwidth D hy D = zt, z be- 
ing the coordination number, e is a reference electronic 
energy and Xj is the j — ih lattice site vector. N is the 
number of molecules in the lattice. It is understood that 
t differs from zero only between first neighbors sites and 
this poses a constraint on the i ^ j sum in the first ad- 
dendum. aj(. ( aic ) creates (destroys) a k-phonon with 
frequency Wk- The lattice dimensionality enters the prob- 
lem through the phonon dispersion relations which have 
been obtained analytically by assuming first neighbors 
pairwise intermolecular forces both along the linear chain 
(ID), the square lattice (2D) and the simple cubic lattice 
(3D) H; 



^ ' M My' 



M M 

= 7^ + 27/3c^ 

Bx,y = 7^(1 + CxCy + SxSy) -I- P^iCx + Cy) 

~^CyCz "t- Sy Sz)) + 2f3'j{cx + Cy + Cz) (0.2) 
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where, Cx = coskx, Cy = cosky, Cz = coskz, = 
sinkx etc. P is the intra-molecular force constant and 
7 is the inter- molecular first neighbors force constant. 
Let's define lu^ = 2(3 /M and ojf = 7/M with M being 
the reduced molecular mass. N is the number of diatomic 
molecules in the lattice and g is the local electron-phonon 
coupling constant. The adiabatic parameter is hCu/D, Q 
being a characteristic phonon frequency which we take 
as the zone center frequency and whose expression is: 
(D^ = Wq + zijj\. 

Throughout this paper we fix Tiloq = IQQnieV and 
t = IbmeV so that the antiadiabatic condition Tiuj / D > 1 
is fulfilled in any dimensionality. Moreover, our pertur- 
bative approach requires the occurence of the condition 
D < g [psil . By applying the Lang-Firsov unitary trans- 
formation pot and the subsequent I/Aq expansion with 
Aq = g'^ / {fibjQD) being the ratio between polaron bind- 
ing energy and electron half bandwidth |30j, H of eq.(l) 
transforms into H = Hq -f Hp with: 



-yy 



exp{{k ■ (r; - r^)) ^ ^ 



k 



1^3 k 



n=0 * k 



(0.3) 



Hq is diagonal except for a second order term in the elec- 
tron density operator which could cause an attractive 
electron-electron interaction The perturbation Hp 
displays the fundamental features of the polaronic quasi- 
particle as the hopping integral narrowing (first factor in 
eq.(3)) plus the peculiar mixing of fermionic and bosonic 
variables. At any electron-phonon interaction vertex m 
(n) phonons can be emitted from (absorbed by) the cloud 
surrounding the electron provided the total crystal mo- 
mentum is conserved. By choosing a transformed ground 
state with no phonons we see that the first order disper- 
sive contribution E^^^ to the ground state polaron band 
arises only from the n — m — term in eq.(3) hence from 
the zero phonon scattering process. In 3D and taking a 
lattice spacing a — jr^ — r, | = 1, one finds 



E^^Hp)— ^'2t{cos px + cos py + cos Pz) 
N 



exp 



Y y i^^kY 



(0.4) 



where the total crystal momentum p coincides with the 
electron momentum due to the absence of self-energy cor- 
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rcctions. The second order perturbative contribution re- 
quires summation over all intermediate states having m 
k-phonons more than the vacuum and one electron on a 
i first neighbor of the j initial site. Moreover, the final 
electronic position / can either coincide with j (this pro- 
cess does not introduce dispersive effects in the polaron 
band) or be a first neighbor of the i site. The latter event 
is clearly dimension dependent: in ID the final site is a 
second neighbor of j, in 2D / can be either a second or 
a third neighbor of j and in 3D, also the fourth neigh- 
bor site can be reached via hopping. While the detailed 
study of these dispersive effects (which can become rele- 
vant in adiabatic conditions) is postponed to a next pa- 
per we turn now to compute the polaron effective mass. 
It should be remarked that the second order corrections 
decrease the mass values calculated in first order pertur- 
bative theory which therefore should be meant to provide 
upper bounds for the polaron mass. 



III. POLARON EFFECTIVE MASSES 

The polaron mass m* can be obtained according to the 
definition 

m* _ zt 

^0 ^ V2i?(p)|p=o ^ ' 

where mo is the bare band mass and the dispersive po- 
laron band is given by eq.(4). The polaron binding energy 
has obviously no p-dispersion. Then, m* /nio is at first 
order independent of t and coincident simply with the re- 
ciprocal of the band narrowing factor. This picture holds 
in the strong coupling regime here assumed. In Fig. 1(a), 
the ratios m*/mo for the ID, 2D and 3D are computed 
versus the first neighbor intcrmolccular force constant 
wi . While the polaron masses strongly depend on the di- 
mensionality d and are very large at small wi, the ratios 
become essentially d independent in the upper portion of 
the parameter range and tend to converge to 2. The value 
of the polaron binding energy Aq > 1 signals that the en- 
ergy gain associated with the lattice deformation is larger 
than the kinetic energy due to the tight binding hopping 
in the bare band. Therefore it is energetically convenient 
to the electron to be dressed by the phonon cloud and be- 
come a quasiparticle. Actually, in antiadiabatic regimes, 
the more restrictive condition concerning the lattice de- 
formation ao = g/ {hu)o) > 1 needs to be fulfilled to guar- 
antee that our quasiparticle is a small polaron. While Aq 
and ao refer to a system with dispersionless phonons it 
is clear that both polaron binding energy and lattice de- 
formation parameter will change after switching on the 
intermolecular interactions. The role of the intermolecu- 
lar couplings is not simply that of increasing the charac- 
teristic phonon frequency but rather that of establishing 
the correct Holstein model dependence on dimensional- 
ity. Ignoring the intermolecular couplings would yield 
the ID polaron band AEid larger than the 2D polaron 
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band AE2D and AE2D > AE^d which is clearly unphys- 
ical since the polaronic wave functions overlap is larger 
in higher dimensionality. This wrong trend would hold 
for any value of the intramolecular frequency loq. Then, 
as observed by Holstein himself in his original paper pT| , 
the phonon dispersion is a vital ingredient of the theory 
and this observation motivates our numerical investiga- 
tion. 

Because of our definitions Aq does not depend on d 
whereas ao is c>c x/d. In Fig. 1(b) we see that the lattice 
deformation a = N^^g^^Tiuj]^^^ is in all dimensions a 
decreasing function of the intermolecular force constant 
and, in ID, the system does not fulfill the small polaron 
condition at the largest wi values. This case has been pre- 
sented to point out how the starting condition aj^ = 1.3 
sets the ID system rather in an intermediate coupling 
regime where a broadening of the polaron size can take 
place [^,^ . Under these conditions the same perturba- 
tive method based on the Lang Firsov transformation be- 
comes questionable . Below tDi — 48mey the polaron 
bandwidth inequalities AE^jj > AE2D > AEuj are not 
satisfied as expected on general grounds hence, the 
dispersionless and the weakly dispersive Holstein Hamil- 
tonian yield erroneous estimates of the effective masses. 
The straight line in Fig. 1(b) marks therefore the lower 
bound for the intermolecular coupling which guarantees 
the validity of the model. 
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FIG. 1. (a) One, two and three dimensional polaron masses 
(in units of the bare band mass) versus the first neighbors in- 
termolecular energy. The dispersionless polaron binding en- 
ergy is 5.3, in units of the bare electron kinetic energy D = zt. 
(b) One, two and three dimensional lattice deformations ver- 
sus the first neighbors intermolecular energy. The are the 
lattice deformation values in a dispersionless model {to-i = 0) . 

= A&meV marks the lower bound for the validity of the 
model (see text). 
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FIG. 2. (a) As in Fig. 1(a) but with a dispersionless polaron 
binding energy Ao — 10.9. (b) As in Fig. 1(b) but with larger 
lattice deformation values. The lower bound for the validity 
of the model is set at a>i = 59meV. 

Increasing the electron-phonon coupling, see Fig. 2(a), 
leads to a strong mass enhancement (particularly in 3D) 
at small oJi while the mass ratios converge to 5 at large 
intermolecular coupling strcnghts. Note (Fig. 2(b)) that 
the polaron is small in all d throughout the whole loi 
range hence the Lang Firsov method works well in this 
case. The threshold value for the validity of the Hol- 
stein model is set here at u^i = 59m.eV. The inequalities 
m^jj > > mljj keep on being satisfied for a por- 

tion of u)i values above the threshold before convergence 
is achieved but second order perturbative terms (being 
larger in higher dimensionality) are expected to correct 
partly this trend. Figs. 3 show that a stronger e-ph cou- 
pling, with Ao = 21.3, yields a mass ratio of ~ 25 and 
shifts the threshold u)i at 65meV pointing out the rela- 
tionship between features of the phonon spectrum and 
strength of the g coupling. Also in this case the po- 
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laron size remains small throughout the whole loi range 
(see Fig. 3(b)) thus confirming the rehabihty of the Lang 
Firsov method in a strong e-ph coupUng regime with an- 
tiadiabatic conditions. 



FIG. 3. (a) As in Fig. 1(a) but with a dispcrsionloss polaron 
binding energy Ao = 21.3. (b) As in Fig.l(b) but with larger 
lattice deformation values. The lower bound for the validity 
of the model is set at a)i = 65meV. 

Next, I have varied g in the range 1 - 4 (in units of 
hojo and found, at any g, the minimum intermolecu- 
lar coupUng Oji{g) at which the bandwidth inequalities 
^EsD > Ai?2Z5 > Ai?i£) are satisfied. This criterion 
yields an empirical relation, (jji{g) — (Z'i(l)(l + ln{g)), 
which allows one to obtain a reliable estimate of the po- 
laron effective mass. In Fig. 4 the ID mass ratio is plot- 
ted versus the dimensionless g/{huJo) both in the first 
and second order of perturbative theory: it turns out 
that second order corrections are negligible in ID sys- 
tems with antiadiabatic conditions as those we have as- 
sumed. I want to point out that the mass values reported 
in Fig. 4 correspond, at any g, to the minimum toi (the 
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threshold) which ensure the smallness of the ground state 
polaron. Then, they are upper bounds for the ID mass 
in the sense that the presence of larger intermolecular 
forces would yield lighter mass values. As expected on 
general grounds |^,|5| the small polaron solution is the 
ground state of the discrete Holstein model in the inter- 
mediate to strong e-p/i-coupling regime here considered 
while, by decreasing the coupling, a continuous cross over 
to large polaron solutions can take place in ID We 
have however seen (Fig. 1(b)) that the dispersive features 
of the phonon spectrum could affect the transition by in- 
ducing a spreading of the lattice deformation. Anyway, 
the smoothness of our m* versus ao curve (persisting 
also in the lower ao range not displayed in Fig. 4) con- 
firms that no self trapping is found in ID antiadiabatic 
regimes whereas recent variational [ p6[ and perturbative 
|33| investigations signalling a rapid growth of m* vs. 
e-ph coupling support the existence of the self trapping 
transition between polaron states of different structure 
in ID adiabatic systems. In any case, phase transitions 
in Holstein models are ruled out being the ground state 
energy an analytic function of the e-ph coupling. 



FIG. 4. One dimensional polaron mass (in units of the 
bare band mass) versus the lattice deformation parameter. 
Both the first and second order perturbative results are dis- 
played. 



IV. CONCLUDING REMARKS 



I have presented the first results of a perturbative ap- 
proach to the polaron problem which focusses on the lat- 
tice dimensionality effects. Having chosen the antiadia- 
batic regime of the Holstein Molecular Crystal Model we 
are confident of the accuracy of the first order perturba- 
tive theory for one dimensional systems with strong e-ph 
coupling whereas some significant second order correc- 
tions can occur in higher dimensionality |^^. While a 
previous work had shown that a dispersionless Hol- 
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stein model leads to i) erroneous estimates of the polaron 
bandwidth versus dimensionaUty and ii) unphysical di- 
vergences in the site jump probabihty the present 
study reveals that the polaron effective mass is in all di- 
mensions very sensitive to the strength of the forces which 
tie the molecules in the lattice. We obtain polaron masses 
between 2 and 25 times the bare band mass by varying 
the e-ph coupling in the range (~ 1 — 2.5) and these val- 
ues become essentially dimension independent when the 
intermolecular forces are sufficiently strong. The molec- 
ular lattice structure has been described by means of a 
single parameter, the first neighbor intermolecular cou- 
pling, being understood that the range of the interactions 
should be extended in real systems if least squares fitting 
of the experimental phonon frequencies can provide ef- 
fective values for the next neighbors and long range force 
constants ^8|] . The antiadiabatic regime with strong e-ph 
coupling ensures the validity of the quasiparticle picture 
for the small polaron nonetheless we have seen that some 
broadening of the phonon cloud can arise at intermedi- 
ate e-ph couplings for strong values of the intermolecular 
forces with consequent lowering of the lattice deforma- 
tion parameter. This interesting feature suggests that 
the intermolecular forces influence the quasiparticle size 
and, incorporating the effects of the e-ph coupling, have 
a role in driving the continuous transition between large 
and small polaron. 
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